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POPULATION MONITORING OF PLAIN PIGEONS IN 
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ABSTRACT.—The Plain Pigeon ( Columba inornata) was listed as endangered in Puerto Rico in 1970. During 
1986-1992 and 1997-2001, we collected point transect survey data to estimate density, population size, and 
rate of change. Density and population size estimates increased between 1986-1992 and 1997-2001. With a 
mean density of 0.25 ± 0.06 SE individuals/ha and a mean population size in the area of the surveys of 3,746 
± 892 SE individuals during 1997-2001, we believe that the status of the Plain Pigeon is not as precarious as 
it was during 1986-1992, when mean density was 0.02 ± 0.003 SE individuals/ha and mean population size in 
the area of the surveys was 218 ± 42 SE individuals. However, Plain Pigeons are not widely distributed and 
the loss and fragmentation of second growth forests combined with the effects of hurricanes and other factors 
may cause their extinction. Because Plain Pigeons have a spatially clumped distribution, we recommend sampling 
at least 1,195 points during peak nesting activity (March through June) throughout the island, with at least 526 
points covering areas of abundance in eastcentral Puerto Rico, to monitor population changes and evaluate the 
effectiveness of management actions. Received 24 June 2002, accepted 11 December 2002. 


Plain Pigeons (Columba inornata) are en¬ 
demic to the Greater Antilles, where they are 
threatened or endangered mainly by habitat 
loss and hunting (Wiley 1985, Perez-Rivera 
1990, Strong and Johnson 2001). Plain Pi¬ 
geons were common and widespread in Puerto 
Rico during the 1800s (Gundlach 1878, Wet- 
more 1916) but became nearly extinct be¬ 
tween the 1920s and 1930s (Danforth 1931). 
They reportedly were hunted between the 
1940s and 1960s, and a small population was 
found at the municipality of Cidra in eastcen¬ 
tral Puerto Rico in 1963 (Leopold 1963). Plain 
Pigeons endured a population bottleneck be¬ 
tween the 1920s and 1960s (Miyamoto et al. 
1994) and were listed as endangered in 1970 
(U.S. Fish and Wildlife Service 1982). 

As with many species listed under the En¬ 
dangered Species Act (Tear et al. 1995), 
guesses and best estimates were used for list¬ 
ing Plain Pigeons. Between 1973 and 1983, 
counts conducted at Cidra fluctuated from 52- 
116 individuals (Wiley 1985), and the popu¬ 
lation was thought to be <200 individuals 
mainly restricted to eastcentral Puerto Rico 
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(Perez-Rivera 1977, 1981). Ad hoc estimates 
suggested that a population increase occurred 
between the 1970s and 1990s. For example, 
Lowe et al. (1990) reported a population of 
about 250 individuals, del Hoyo et al. (1997) 
reported about 300 individuals, Birdlife Inter¬ 
national (2000) reported about 700 individu¬ 
als, and Oberle (2000) reported about 2000 
individuals. 

Reliable population parameter estimates are 
needed to establish recovery goals (Tear et al. 
1995) and to assess population viability under 
adverse deterministic and stochastic process¬ 
es, such as habitat loss and catastrophic 
weather (Lande 1993, Saether et al. 1998). 
Standard sample surveys were conducted in 
Puerto Rico during 1986-1992 and 1997- 
2001 to estimate density, population size, and 
rate of change over time (Rivera-Milan 1992, 
1993, 1995a; Martinez et al. 2001). Based on 
these parameter estimates, we assessed the 
population status of Plain Pigeons and make 
recommendations to continue monitoring 
trends as part of recovery efforts (U.S. Fish 
and Wildlife Service 1982, Puerto Rico Dept, 
of Natural and Environmental Resources 
2000). 

METHODS 

Study area .—Puerto Rico, the smallest 
(880,190 ha) and easternmost (18° 15' N, 66° 
30' W) island of the Greater Antilles, lies 
within the subtropical region, and has six life 
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zones (Ewel and Whitmore 1973). Its forest 
cover increased from about 6% during the 
1940s (Koenig 1953, Birdsey and Weaver 
1982) to about 42% during 1991-1992 (Hel- 
mer et al. 2002). About 5% of the island area 
is protected, but rapid and widespread devel¬ 
opment is causing the loss and fragmentation 
of second growth forests (Birdsey and Weaver 
1982, 1987; Ramos-Gonzalez 2001; Helmeret 
al. 2002). 

To collect survey data, we stratified the is¬ 
land according to the area of its three major 
life zones (dry zone = 121,640 ha, moist zone 
= 532,610 ha, and wet zone = 212,480 ha; 
Ewel and Whitmore 1973). Most of the forests 
in these life zones are characterized by second 
growth vegetation at different stages of suc¬ 
cession (Birdsey and Weaver 1982). For ex¬ 
ample, early succession second growth forests 
in the moist zone are characterized by species 
such as trumpet tree ( Cecropia schreberiana ), 
American muskwood ( Gucirea guidonea), Af¬ 
rican tuliptree ( Spathodea campamdata ), and 
guaba (Inga vera). 

Point transect surveys. —During each year 
of the study we conducted point transect sur¬ 
veys during peak nesting activity (March 
through June) on 8-km routes along secondary 
(paved) and tertiary (unpaved) roads covering 
all habitats used by columbids in the life 
zones. Depending upon topographic features 
(e.g., winding versus straight roads), the 
routes had 6-11 points (stations), with the first 
point of each route placed randomly and the 
others systematically at 800-m or 1,600-m in¬ 
tervals to minimize the chances of counting 
the same individuals more than once. We 
marked point centers so they could be found 
easily at all times. We used range finders to 
measure distances from point centers to Plain 
Pigeons detected singly or in flocks (clusters 
>2 individuals). Because we recorded dis¬ 
tances to all columbids detected from point 
centers, distances were grouped into catego¬ 
ries (0-15, 15-30, 30-45, 45-60, 60-90, 90- 
120, 120-180, 180-240, 240-340. 340-440, 
and >440 m) to secure data quality (Buckland 
et al. 1993). 

We estimated density using the formula: 

D = nh(0)s/27rk, 

where D = the number of individuals/ha; n = 
the number of individuals counted singly or 


in clusters; h(0) = the slope of the probability 
density function of detection distances, esti¬ 
mated at zero distance; s = mean cluster size; 
and k = the number of points sampled. At 
each point, two observers counted for 6 min 
all the Plain Pigeons seen or heard. We re¬ 
corded distances to individuals flying over¬ 
head only if they were seen landing and their 
locations could be fixed before they started 
moving again. We recorded the distance from 
the center of a point to the center of a cluster 
of individuals or single individual up to 440 
m. Detections beyond 440 m were recorded 
but not used for density estimation. Because 
columbids nest year-round (Wiley 1991; Ri- 
vera-Milan 1996, 2001; Rivera-Milan et al. 
2003) and share a typical pattern of nest at¬ 
tendance (with males on the nest from mid¬ 
morning to midafternoon or early evening, 
and females on the nest from late afternoon 
until midmorning of the next day; Blockstein 
and Westmoreland 1993), we conducted the 
surveys from 06:00-10:30 AST to include 
breeding and nonbreeding individuals of both 
sexes. 

We used program DISTANCE (Thomas et 
al. 2001) to analyze the point transect survey 
data. We evaluated the fit of detection models 
(uniform, half-normal, and hazard-rate key 
functions with cosine, simple polynomial, and 
Hermite polynomial parameter adjustments) 
to data with x 2 goodness-of-fit tests. Model 
selection was based on Akaike’s Information 
Criteria (AIC or AIC c ; Burnham and Ander¬ 
son 1998). Generally, data beyond 180 m 
were truncated to remove outliers and im¬ 
prove model fitting. If the regression of nat¬ 
ural logarithm of cluster size (ln[sj) and de¬ 
tection distance (r) was significant at a = 
0.15, we used the size-bias regression method 
to estimate mean cluster size. Encounter rate 
variance was estimated empirically. Based on 
data collected in eastcentral Puerto Rico from 
1997-2001 (k = 250) and throughout the is¬ 
land from 1998-2001 (k = 765-1,065), we 
calculated the number of points (p,) needed to 
obtain a coefficient of variation (CV) of 20% 
for estimated density (Buckland et al. 1993: 
304-307). 

We used a one-tailed z- test (Thompson et 
al. 1998) to test the hypothesis that log-trans¬ 
formed estimates of density and population 
size of Plain Pigeons in eastcentral Puerto 
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1986-1992 1997-2001 



FIG. 1. Box plots of (A) density (D = the number 
of individuals/ha) and (B) population size (N = the 
number of individuals in the area of the surveys) es¬ 
timates based on point transect survey data collected 
for Plain Pigeons in eastcentral Puerto Rico from 
March through June, 1986-1992 (k = 187) and 1997- 
2001 {k = 250). Percentile ranks (10th, 25th, 50th, 
75th, and 90th) are shown, with the top and bottom of 
each box plot representing the 50th percentile, the line 
in the middle representing the median, the line extend¬ 
ing above representing the 90th percentile and the line 
extending below representing the 10th percentile, and 
the circles extending from the lines representing values 
below the 10th percentile and above the 90th percen¬ 
tile. 


Rico were lower during 1986-1992 (k = 187) 
than during 1997-2001 (k = 250). To test for 
positive trends ((3, > 1) in log-transformed es¬ 
timates, we used simple linear regression (Nur 
et al. 1999). Residual plots and residual sta¬ 
tistics indicated that the assumptions of re¬ 
gression were met [e.g., ln(population size): 
Durbin-Watson d = 1.6, P > 0.05; serial au¬ 
tocorrelation = 0.2, Box Ljung Q = 0.7, P = 
0.40]. 

RESULTS 

Plain Pigeon density estimates ranged from 
0.01-0.03 individuals/ha (x = 0.02 ± 0.003 
SE) and population size estimates in the area 
of the surveys ( krrr 2 = 1 1,374 ha) ranged from 
60-377 individuals (x = 218 ± 42 SE) in 
eastcentral Puerto Rico during 1986-1992. 
Density estimates ranged from 0.10-0.40 in¬ 
dividuals/ha (x = 0.25 ± 0.06 SE) and pop¬ 



FIG. 2. Distance (r) and detection probability 
(g[>]) of Plain Pigeons based on point transect survey 
data collected on Puerto Rico during April through 
June, 2001 (n = 64 and k = 1,065). The half-normal 
key function + 2-term cosine parameter adjustments 
provided the most parsimonious fit to the data. 


ulation size estimates in the area of the sur¬ 
veys (kirr 2 = 15,205 ha) ranged from 1,497— 
6,106 individuals (x = 3,746 ± 892 SE) dur¬ 
ing 1997-2001. Estimated density ( z = —2.5, 
P = 0.007) and population size (z = —2.8, P 
= 0.002) were lower during 1986-1992 than 
during 1997-2001, respectively (Fig. 1). Rate 
of change was positive and significant (den¬ 
sity: p, = 0.29 ± 0.07^ SE, /, 2 = 3.9, P = 
0.001; population size: |3, = 0.33 ± 0.09 SE, 
t ]2 = 3.7, P = 0.002). 

Plain Pigeons rarely were detected outside 
24 routes (k = 187-250) in eastcentral Puerto 
Rico. In general, detection probability re¬ 
mained high in the first five distance catego¬ 
ries (effective detection distance: x — 88 m ± 
5 SE), cluster size (x = 1.8 ± 0.1 SE) was 
not related to detection distance, and the half¬ 
normal key function with cosine parameter 
adjustments fitted the data. For example, 64 
detections were made in 1,065 points sampled 
from April through June, 2001, which resulted 
in a mean effective detection distance of 59 
m ± 5 SE, a mean cluster size of 1.03 ± 0.02 
SE (regression of hi[s] against r: B, = 0.0003 
± 0.0003 SE, t 62 = 0.9, P = 0.83), and a half¬ 
normal key function with cosine parameter 
adjustments (m = 2) best fitting the data (AIC c 
= 245.7; x 2 4 = 2.3, P = 0.67; Fig. 2). 

During 1998-2001 (k = 765-1,065), the 
coefficient of variation of density estimates 
ranged from 25-50% (x = 35% ± 0.5 SE) and 
the dispersion parameter estimate ( b ) ranged 
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from 1.9-10.7 (T = 6.9 ± 1.8 SE). Because 
the spatial distribution of Plain Pigeons was 
clumped (b > 1), we needed to sample at least 
1,195 points throughout the island, with at 
least 526 points in eastcentral Puerto Rico, to 
obtain a desired coefficient of variation of 
20% for estimated density. 

DISCUSSION 

In Puerto Rico, Plain Pigeons are multiple- 
brooded habitat generalists that feed on the 
fruits of a wide variety of plant species, main¬ 
ly in mesic second growth forests (Perez-Ri- 
vera 1978). Because density and food abun¬ 
dance are positively related (Rivera-Milan 
1992, 1996, 2001; Bancroft et al. 2000; Ri¬ 
vera-Milan et al. 2003), we suggest that the 
recovery of second growth forests, particular¬ 
ly during the 1970s and 1990s, increased the 
quantity and quality of foraging and nesting 
habitats, which in turn elicited a positive pop¬ 
ulation response in eastcentral Puerto Rico be¬ 
tween 1986-1992 and 1997-2001. Forest area 
has increased since the 1940s because mar¬ 
ginally productive pasture and cropland were 
abandoned as Puerto Rico gradually became 
industrialized (Koenig 1953, Birdsey and 
Weaver 1982, Helmer et al. 2002). 

With a mean density of 0.25 individuals/ha 
and a mean population size in the area of the 
surveys of 3,746 individuals during 1997- 
2001, we believe that the status of the Plain 
Pigeon is not as precarious as it was during 
1986-1992, when mean density was 0.02 in¬ 
dividuals/ha and mean population size in the 
area of the surveys was only 218 individuals. 
Our optimism is tempered, however, because 
Plain Pigeons are not widely distributed and 
unmitigated development is causing major 
land cover changes (Ramos-Gonzalez 2001, 
Helmer et al. 2002), which may be affecting 
the reproduction of Plain Pigeons through the 
loss and fragmentation of second growth for¬ 
ests in eastcentral Puerto Rico (Wiley 1985, 
Perez-Rivera 1990, Rivera-Milan et al. 2003). 
Habitat changes may interact with catastroph¬ 
ic weather, food availability, nest predation, 
illegal hunting, and poaching of nestlings to 
reduce persistence, even when the population 
can be increasing exponentially (Mangel and 
Tier 1994). Illegal hunting, poaching of nest¬ 
lings, and other forms of human-induced dis¬ 
turbance still occur (Ruiz-Lebron et al. 1995) 


but not as frequently as during the 1970s (U.S. 
Fish and Wildlife Service 1982, Wiley 1985). 

Hurricanes can affect plant and animal pop¬ 
ulations through multiple effects (Wiley and 
Wunderle 1994, Wunderle 1995, Lomascolo 
and Aide 2001). During 1986-2001, two hur¬ 
ricanes struck Puerto Rico: Hugo (a Category 
5 hurricane based on the Saffir-Simpson scale) 
in September 1989 and Georges (Category 3) 
in September 1998. Hugo had stronger winds 
but Georges affected a larger portion of the 
island (for more information, see http:// 
www.geocities.com/huracanado 1/). Estimates 
of density and population size in November 
1998 (D = 0.01 ± 0.01 SE and N = 310 ± 
182 SE) were as low as during March through 
June, 1986-1992, remained depressed during 
February through October, 1999 (D < 0.10 ± 
0.03 SE and N < 1,497 ± 523 SE), and re¬ 
bounded to pre-hurricane levels during March 
through June, 2000-2001 (D > 0.12 ± 0.04 
SE and N > 1,772 ± 644 SE). 

As with other columbids (Rivera-Milan 
1995b, Wunderle 1995), Plain Pigeons 
showed resiliency (Pimm 1991), and through 
successful reproduction responded rapidly to 
forest regeneration and increased food avail¬ 
ability (Rivera-Milan 1996, 2001; Rivera-Mil¬ 
an et al. 2003). However, a hurricane with 
Hugo’s wind force striking the island as did 
Georges can be devastating for Plain Pigeons. 
We hypothesize that extensive deforestation 
from agriculture (Koenig 1953, Birdsey and 
Weaver 1982), unregulated hunting pressure 
(Wiley 1985), and the passing through most 
of the island of hurricanes San Felipe (Cate¬ 
gory 5) in 1928 and San Ciprian (Category 3) 
in 1932, interacted to endanger an already re¬ 
duced and patchily distributed population. 

Reliable population monitoring data are 
needed to evaluate the effectiveness of man¬ 
agement actions (Gibbs et al. 1999) and de¬ 
termine if recovery goals have been reached 
(Tear et al. 1995). Distance sampling theory 
provided an adequate framework to survey 
columbids using roads. In the case of the Plain 
Pigeon, this framework needs to reflect ad¬ 
justments resulting from their spatially 
clumped distribution. This distribution may be 
the result of social behavior, reproductive ad¬ 
aptations, and habitat use with respect to the 
spatial and proportional distribution of forag¬ 
ing and nesting habitats. Plain Pigeons can 
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disperse long distances in search of food 
(Ruiz-Lebron et al. 1995), but nesting has not 
been reported outside eastcentral Puerto Rico 
(Rivera-Milan 1996, 2001; Rivera-Milan et al. 
2003). Because second growth forests have 
similar vegetation structure and composition 
in the moist zone, we do not know why they 
are not more widely distributed. 

A spatially clumped distribution means that 
many points will have no detections whereas 
a few points will have detections of individ¬ 
uals or clusters of individuals. Depending 
upon the degree of clumping, sampling sev¬ 
eral hundred or thousand points may be re¬ 
quired to obtain reasonably precise density es¬ 
timates [CV(D) < 20%]. Currently, 1,065 
points are sampled on Puerto Rico annually. 
This sampling effort is enough to obtain pre¬ 
cise density estimates for abundant and widely 
distributed columbids, such as Scaly-naped 
Pigeons (Columba squamosa ), Zenaida Doves 
(Zenaida aurita), and White-winged Doves 
(Z. asiatica). So, given the cost of conducting 
wide-ranging surveys, we suggest sampling at 
least 1,195 points throughout the island, with 
at least 526 points covering eastcentral Puerto 
Rico. Surveys should be conducted from 
March through June to cover the peak of nest 
density of Plain Pigeons and other columbids 
(Rivera-Milan 1996, 2001; Rivera-Milan et al. 
2003). Surveys conducted from March 
through June (CV: x = 33% ± 2 SE) were 
less variable than surveys conducted from 
July through February (CV: x = 47% ± 7 SE) 
because of changes in flocking behavior (Ri¬ 
vera-Milan 1992, 1993, 1995a; FFR-M un- 
publ. data). 

Plain Pigeons were detected easily and 
grouped into distance categories ranging from 
0-15 m to 340-440 m. In general, detection 
probability remained high near point centers 
and decreased with distance, cluster size was 
not a function of detection distance, and de¬ 
tectability followed a half-normal process 
(Buckland et al. 1993). Plain Pigeons were de¬ 
tected near (0-30 m) and far (>120 m) from 
roads because they frequently perched at can¬ 
opy level and moved conspicuously within 
and between forest fragments. To be conser¬ 
vative, however, we did not estimate popula¬ 
tion size beyond the surveyed area ( kirr 2 ) be¬ 
cause placement of sampling units on or near 
roads may not be representative of their abun¬ 


dance in off-road habitats (Buckland et al. 
1993, Thompson 2002). Further research is 
needed to validate sample inferences from on¬ 
road to off-road habitats (Keller and Scallan 
1999). 

Islandwide surveys would be impossible 
without the use of roads. Secondary and ter¬ 
tiary roads provide access to all habitat types 
used by columbids, ranging from upper mon¬ 
tane wet forests to coastal dry forests. In 
Puerto Rico, Plain Pigeons behave as edge 
species (Yahner 1988), nesting, foraging, and 
roosting in trees at or near roads. Thus, we 
believe that point transect surveys along roads 
provided adequate coverage of their habitats 
and daily activities, and we recommend com¬ 
bining distance sampling with telemetry and 
geographic information systems to link de¬ 
mographic and spatial data. These data should 
be used to establish specific population and 
habitat recovery goals (Tear et al. 1995) for 
the downlisting and delisting of Plain Pigeons 
in Puerto Rico. 
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